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. Abstract. The supranova model for gamma-ray bursts, originally pro- 

CN I posed by Vietri &: Stella (1998), has several unique features that make 

• it attractive for the interpretation of GRBs and their afterglows, and it 

D I has emerged as a promising candidate for modeling the evolution and 

environment of at least some GRB sources. This review summarizes the 
^ I model and its key observational implications, assesses its strengths and 

potential weaknesses, and outlines paths for future observational and the- 
oretical work. 
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' 1. Introduction 

cn ' 

I There is mounting evidence that the progenitors of long-duration 2 s) gamma- 

. ray bursts (GRBs) are massive stars that collapse to form stellar-mass black holes 

Qh| or strongly magnetized and rapidly rotating neutron stars. Up until recently, 

Q ■ the prevalent scenario has been the "collapsar" model, in which the GRB event 

j-l I directly follows the stellar collapse (e.g., Woosley, these Proceedings). This pic- 

t/^ ■ ture is consistent with the identification in a few GRBs of a possible signature 

of a nearly simultaneous supernova Ib/Ic explosion. In one commonly discussed 
realization of this scenario, the duration of the burst is determined by the accre- 
I tion time of the unexpelled stellar material through the debris disk that forms 

^ • around the newly created central black hole (BH). 

Vietri & Stella (1998) proposed an alternative to the collapsar picture, 
dubbed the "supranova" model, in which the supernova (SN) explosion initially 
results in the formation of a comparatively massive, magnetized neutron star 
(NS) endowed with rapid rotation. This "supramassive" NS (SMNS) is envi- 
sioned to gradually lose rotational support through a pulsar-type wind^ until it 
eventually becomes unstable to gravitational collapse (leading to the formation 
of a BH and the triggering of a GRB). The original motivation for the supranova 
model was the desire to identify a comparatively "baryon clean" environment 
in which the high (^ 100) Lorentz factors inferred in the prompt high-energy 
emission region could be attained. (In the proposed picture, the gas around the 
GRB progenitor is swept out by the expanding SN ejecta over the spindown time 
of the SMNS.) As discussed below, subsequent considerations of this scenario 



^It is assumed in what follows that angular momentum loss by gravitational waves remains 
relatively unimportant. 
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have revealed that it also has the potential of providing a natural interpretation 
of several key observational features of GRB afterglows. 

2. The Supranova Scenario 

2.1. SMNS Evolution 

An SMNS is a gcncral-rclativistic equilibrium configuration of a uniformly ro- 
tating NS whose mass exceeds (by up to 20%, for realistic equations of state) 
the maximum mass of a nonrotating NS (e.g., Cook et al. 1994; Salgado et al. 
1994).^ An SMNS that loses angular momentum and energy adiabatically while 
conserving its total baryon mass follows an evolutionary sequence that brings 
it to a point where it becomes unstable to axisymmetric (and possibly first to 
nonaxisymmetric) perturbations, which result in collapse to a BH. 

For an SMNS of mass M^,, radius i?*, angular velocity 51*, disposable (before 
the SMNS collapses) rotational energy AE,ot = aGMli\/2c k. 10^^ gj.gs, and 
surface magnetic field 5*, the spindown time due to a pulsar-type wind is 

^\ / M* \2 / \-6 / \ -3 / 

O.sj UmoJ Uskmj Uo^s-ij UOI^G; 

which can range from several weeks to several years for plausible values of M*, 
i?*, and (~ 10^2 _ iQi3 Q) 

One possible scenario of how the formation of a BH could trigger the GRB 
outburst is that a transient debris disk comprising the outer layers of the SMNS 
is left behind when the bulk of the SMNS collapses, and that the stellar magnetic 
field threading this material is amplified by differential rotation to 10^^ G and 
drives a Poynting flux-dominated outflow from the disk surfaces (e.g., Meszaros 
& Rees 1997; Vlahakis &; Konigl, these Proceedings). 

2.2. Theoretical Foundations of the Model 

Only a small fraction of newly formed neutron stars could belong to the supra- 
massive class: this is consistent with the fact that GRBs are rare events, but so 
far our knowledge of how a rotating star collapses is not sufficient to make an 
a priori determination of this fraction. However, the results already available 
suggest that this fraction may well be nonnegligible. For example, numerical 
simulations of nonrotating stars have been used to infer that the distribution of 
NS masses produced through core collapse may be flat beyond the 1.2 — 1.6 Mq 
range where it peaks (Fryer & Kalogera 2001). Not unexpectedly, when rota- 
tion is included, the simulations yield more massive remnants (Fryer & Heger 



^Differential rotation may support a significantly larger mass against gravitational collapse, 
allowing the formation of a hypermassive NS (e.g., Baumgarte et al. 2000). However, a dif- 
ferentially rotating NS is estimated to evolve to a uniform rotation profile (through magnetic 
field amplification and viscosity) on a timoscalc much shorter than the spindown time of a 
uniformly rotating star (e.g., Shapiro 2000). Therefore, in contrast with supramassive config- 
urations, hypermassive neutron stars above the supramassive mass limit are not expected to 
exhibit significant delays between NS formation and collapse to a BH. 
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2000). The NS angular velocity distribution derived from stellar-evolution cal- 
culations is sensitive to the adopted assumptions. These studies indicate that, 
in the absence of magnetic angular-momentum transport, newly formed neutron 
stars will likely rotate near breakup (Heger et al. 2000). Recently proposed pre- 
scriptions for the inclusion of magnetic effects (e.g., Spruit 2002) imply strong 
angular-momentum transport and comparatively low initial values of fi* (Heger 
& Woosley 2002), but their validity remains uncertain. In fact, scaling argu- 
ments from white-dwarf systems suggest that corc-to-cnvclopc angular momen- 
tum transport in pulsar progenitors may be rather inefficient (Livio & Pringle 
1998). 

Another relevant question is whether a debris disk — perhaps the most 
likely origin of a GRB outflow - is left behind when the SMNS eventually col- 
lapses to a BH. Although a collapsing supramassive compact object with a soft 
equation of state could leave behind up to 10% of its mass in the form of 
a rotationally supported disk (e.g., Shapiro &; Shibata 2002), realistic SMNS 
equations of state might be too stiff in this regard, with the implication that the 
entire SMNS would collapse (e.g., Shibata ct al. 2000). However, comparatively 
stiff SMNS configurations could reach the mass-shedding limit before they be- 
come dynamically unstable (e.g., Cook et al. 1994), so in principle such stars 
might also give rise to a BH surrounded by a highly magnetized debris disk.^ 

3. GRB Afterglows in Pulsar-Wind Bubbles 

The supranova scenario also has important implications to GRB afterglows. In 
particular, the bubble created by the SMNS pulsar-type wind provides an envi- 
ronment that can account for both the high fraction (eg) of the internal energy 
residing in rclativistic electrons and positrons and the high magnetic-to-internal 
energy ratio (es) that have been inferred in the afterglow emission regions of a 
number of sources through spectral modeling. The deduced values of Cg (^ 0.1; 
e.g., Panaitescu & Kumar 2002) are hard to explain in a standard [interstellar- 
medium (ISM) or stellar-wind] environment, especially for long-lasting after- 
glows. The values of are typically inferred to be not much smaller than eg 
(for example, spectral fits to the GRB 970508 afterglow have yielded esti- 
mates in the range 0.01 — 1). These estimates are again problematic for GRB 
outflows that propagate into a standard ambient medium, since €b will only 
be ~ lO"^'^ for a Milky Way-like ISM and at most ^ 10~^ for a magnetized 
stellar- wind environment. 

The comparatively high values of eg and €b can, however, arise naturally 
in a pulsar- wind bubble (PWB), since the composition of rclativistic pulsar 
winds is likely dominated by an e+e~ component, and since such winds are 
often characterized by a high magnetization parameter aw = B^/AttPwC^ (the 
Poynting-to-particle energy flux ratio, where By, is the magnetic field amplitude 
and pw is the rest-mass density, both measured in the wind frame). This has 
been explicitly demonstrated by Konigl & Granot (2002; hereafter KG02). 



■^Although additional material (e.g., some of the fallback from the original SN explosion) could 
surround the newly formed BH, this gas would not be as highly magnetized as the SMNS 
matter and therefore might not contribute significantly to the GRB outflow. 
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KG02 constructed a PWB model based on a number of simplifications — in- 
cluding spherical symmetry, negligible ambient mass, a steady state, pure e+e^ 
composition, a monoenergctic pair distribution at each radius, and a dominant 
azimuthal magnetic-field component. In their representation, the SMNS-driven 
pulsar-type wind shocks at a radius Rs and fills the volume between Rs and 
the radius Rb of the swcpt-up SN-ejecta shell with a relativistic fiuid consist- 
ing of "hot" e^e~ pairs and electromagnetic fields. The ejecta shell (of mass 
Mcj ~ IOMq) is accelerated by the bubble (becoming compressed and frag- 
mented in the process) and expands at a speed Vb ~ 

The afterglow emission is attributed to the bow shock of a GRB jet that prop- 
agates between Rg and Rb- KG02's choice of wind parameters was guided by 
recent observations of the Crab and Vela SN remnants (SNRs), and their model 
incorporates synchrotron-radiation cooling (parameterized by ai ^ synchrotron 
cooling time in units of Rb/c)'^ and a possible departure from ideal MHD (asso- 
ciated with an "equipartition" upper bound on the electromagnetic-to-thermal 
pressure ratio) in the shocked-wind zone. 

KG02 showed that their model can account for the inferred values of eg and 
eg as well as for the deduced range of ambient densities in GRB afterglows (e.g., 
Panaitescu 8z Kumar 2002). In this interpretation, the ambient density is not 
dominated by the baryon rest mass as in the standard picture, but rather by the 
relativistic inertia of the shocked wind and compressed electromagnetic fields, 
so the effective preshock hydrogen number density is given by 

_ Ap + {B + Ef/Air 

'^H,equiv — 9 ) 

mpC^ 

where p, B, and E are, respectively, the particle pressure and the magnetic and 
electric field amplitudes in the PWB (all measured in the fluid rest frame). By 
using the parameterization raH,equiv(^) r~^, it was found that k takes on a 
range of values (determined by cr^, ai, and the location within the bubble) that 
encompass both a uniform-ISM and a stellar-wind radial density profiles. The 
derived solutions have demonstrated that in the context of this scenario it is 
possible to explain how a GRB with a massive stellar progenitor can produce 
an afterglow that shows no evidence of cither a stellar wind or a high-density 
environment, which has been hard to understand in the standard picture. Ad- 
ditional observational implications of the PWB interpretation are discussed by 
Guetta k. Granot in these Proceedings. 

4. Interpretation of X-ray Features 

As reviewed in one of this Workshop's sessions, several afterglow sources have 
shown evidence for X-ray emission (and absorption) features, which have been 
identified as either Fc or primarily "light" -clement (Mg, Si, S, Ar, Ca) transi- 
tions. Although it has been recognized early on that these features are consis- 
tent with an origin in a large-scale SNR shell and could thus lend support to 



^As was pointed out by Guctta & Granot (2003; hereafter GG03), the effect of cooling on the 
PWB structure is less important if ions dominate the wind particle energy flux. 
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the supranova model, alternative explanations that might be compatible with 
the collapsar model have also been proposed (e.g., Kallman et al. 2001). In 
the following paragraphs I first summarize the interpretation of the reported 
emission features that was outlined in KG02, and then I discuss its implications 

to the supranova scenario. 

In the supranova picture, the emission is induced by continuum irradiation 
from the central region around the time of the burst. Light-travel effects limit 
the solid angle from within which emission is received at time t to A$7(t)/47r = 
[1 — cos6{t)]/2 = ct/2{l + z)Rb, where z is the source rcdshift. For example, in 
the case of GRB 991216 (observed between 37 and 40.4 hr; Piro ct al. 2000), the 
model implies that the detected emission originated at polar angles 9 ~ 27° —28° 
with respect to the jet axis. The emission feature is interpreted as Fe XXV Hea 
from a photoionized gas of ionization parameter ^ = Li/nR^ 10^, where 
is the ionizing continuum luminosity. For Lj = ^ttD'^Fx = 6.1 x 10^^ ergs s~^, 
Mpe ~ 0.1 Mq, and T ss 10^ K, one infers iZ^ « 2 x 10^^ cm. Since the expansion 
time to Rb ~ 10^^ cm at a speed Vb ~ 0.1c is ~ 10^ days, this estimate is 
consistent with the bulk of the SN-ejected radioactive ^^Ni decaying into ^^Fe. 
The estimated density is s» 2 x lO^'^ cm~^, which likely corresponds to clumps 
in the Rayleigh- Taylor unstable SNR shell. If the clumps have a covering factor 
~ 1 and an iron abundance a few times solar, then the implied Thomson optical 
depth tt is ^ 1 and the Fc photoionization optical depth is a few, which are 
optimal for producing high -equivalent- width iron lines through reflection.^ 

The identification of X-ray features in a source like GRB 011211 with "light" 
elements (Reeves et al. 2002) can be interpreted in this picture in terms of a 
bubble characterized by igd ^ weeks (rather than months), in which iron has 
not yet formed in the shell (Granot & Guetta 2002). 

The value of Rb inferred from modeling the X-ray emission feature in GRB 
991216 is smaller (by about an order of magnitude) than the distance deduced 
from spectral fits to the continuum afterglow emission. As was pointed out by 
KG02, a plausible resolution of this apparent discrepancy is that, rather than 
being spherical, the PWB is, in fact, elongated along the SMNS rotation axis. 
(Further discussion of this possibility, as well as additional arguments against 
a simple spherical model, are given in GG03.) Such a morphology is a likely 
outcome of the PWB structure and of its origin. In particular, 

• The supernova explosion of a rapidly rotating progenitor star is likely to be 
nonspherical, and so is also the pulsar wind; 

• Both the SN ejecta and the environment into which the SNR expands (which 
was possibly shaped by earlier episodes of stellar mass loss during the red- 
supergiant and blue-supergiant evolutionary phases of the progenitor star) could 
have a highly anisotropic mass distribution; 

• Magnetic hoop stresses may act to strongly collimate the PWB (and possibly 
also the preceding stellar outflows), an effect enhanced by cooling (e.g.. Chevalier 
& Luo 1994; Gardiner & Frank 2001). 

In this interpretation, a narrow GRB jet moves within the PWB along 
the rotation axis, and the afterglow shock reaches much larger distances from 



^Reflection from a medium with tt S> 1 is, however, also a possibihty; e.g., Ballantyne et al. 
(2002). 
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the center than the spherical portion of the SNR shell from which the X-ray 
emission features originate. This picture could potentially be corroborated by 
observations.^ For example, it implies (on account of light-travel effects) a min- 
imum onset time for the detectability of X-ray emission features. There have 
already been several reported null detections of such features in sensitive X-ray 
observations of GRB sources: the above picture suggests that continuing the 
search in sources like these at later times might be worth while. Another impli- 
cation of the above interpretation is that the polar region of the bubble expands 
faster (by a factor of up to ^ 10) than the lower-latitude SNR shell where the 
emission features are produced. If a dense SN-ejecta clump located near the ma- 
jor axis of the PWB was carried out by the expanding bubble and happened to 
intercept our line of sight, then an absorption feature with an implied blueshift 
significantly larger than that of the emission feature might be detected.^ 



5. Assessment of the Model 

Several recent papers that scrutinized the supranova scenario have drawn atten- 
tion to a number of potential weaknesses of this model. I now briefly summarize, 
and comment on, some of the issues that have been raised. 

• Occurrence of long bursts. It has been argued (e.g., Bottcher & Fryer 2001) 
that the supranova scenario is most likely to produce short bursts, rather than 
the long bursts to which this model has been applied. There are, however, 
potential alternatives to ongoing mass feeding (as in the collapsar model) for 
producing long disk lifetimes. They include a low disk viscosity (e.g., Popham 
et al. 1999; Ruffcrt &: Janka 1999) and a magnetically mediated spin-up torque 
exerted by the BH (van Putten &; Ostriker 2001). 

• Apparent evidence (from the detection of a reddened bump in the GRB 
lightcurve in a few sources) for a SN/GRB coincidence. If this interpretation 

is correct, it argues against a time delay between the SN explosion and the 
GRB event. However, alternative interpretations of these bumps have been pro- 
posed, including illumination of ambient dust by the GRB source (e.g., Esin &i 
Blandford 2000; Waxman & Draine 2000) and (in the context of the supranova 
scenario) thermal emission from a PWB-heated SNR shell (Dermer 2002). Fur- 
thermore, as discussed by GG03, the SN/GRB time delay may span a range of 
values, with a near coincidence (as in the collapsar model) possibly characteriz- 
ing a subset of sources. 

• Disruption of SNR shell. In contrast to the PWB model described in §3, Inoue 

et al. (2003) explored a scenario wherein the pulsar wind disrupts and penetrates 
through the SN-ejecta shell, producing a bubble whose size is about an order 
of magnitude larger than the SNR radius. One can argue, however, that since 
Vb as calculated in the KG02 picture is much less than the wind speed (~ c), 



Given, however, that a variety of factors could in principle prevent the effects discussed here 
from being observed, a failure to detect them in any given source would not necessarily con- 
stitute evidence against this interpretation. 

preliminary report on a tentative detection of such an absorption feature in GRB 011211 
w£is given in F. Frontera's presentation in this Workshop. 
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it is likely that the wind will shock and remain confined within the SNR shell. 
Existing numerical simulations of "plerionic" SNRs (e.g., Jun 1998) seem to be 
consistent with this conclusion, but additional studies that explicitly consider 
the case where the total wind energy greatly exceeds that of the SN and that 
explore strong departures from spherical symmetry would be useful. 

• Apparent need for fine-tuning. McLaughlin et al. (2002) observed that the 
typical shell size 10^^ cm) inferred in the supranova interpretation of the 
detected X-ray emission features translates into a characteristic value of the 
SN/GRB time delay: they questioned why this value would be constant from 
source to source and wondered whether it would be long enough to account for 
the decay of nickel in the SN ejccta into iron. However, as was noted in §4, there 
is now evidence for X-ray features arising in "light" elements rather than iron, 
which can be naturally interpreted in the supranova picture in terms of lower 
values of tg^- As was also noted in §4, the time delays required for iron to become 
dominant are consistent (in view of the ~ 0.1c shell expansion speed estimated in 
the spherical model) with Rg being ^ 10^^ cm. Furthermore, as was pointed out 
by KG02, shells with sizes much larger than 10^^ cm would have densities that 
are too low to give rise to observable features. These considerations can explain 
why the inferred radii of iron-emitting shells are roughly the same. Another fine- 
tuning issue was raised by Woosley et al (2002), who called attention to the fact 
that SMNS configurations correspond to a limited NS mass range, especially in 
comparison with hypermassive neutron stars (in which significant SN/GRB time 
delays are not expected; see footnote 2). Although this statement is factually 
correct, it seems to be more of a description of a constraint that the model needs 
to satisfy (which may be part of the reason why GRBs are rare events) than of 
a serious flaw. 

In summary, even though the ultimate verdict on this model will be provided 
by additional observations and theoretical calculations, it appears that none of 
the objections raised against it so far is compelling. 

6. Conclusions 

The supranova model shows great promise in that it naturally addresses several 
key modeling issues of GRBs and their afterglows: 

• a relatively baryon-free near-source environment (the pulsar-wind cavity) in 
which a high-Lorentz-factor outflow can be accelerated; 

• a pair-rich and potentially strongly magnetized environment further out (the 
pulsar- wind bubble) in which a high-eg and es afterglow emission can arise; 

• a heavy element-enriched, extended dense medium (the clumpy SNR shell) in 
which the X-ray spectral features can be produced. 

The model can readily account for a variety of observed behaviors in GRB 
afterglows, which may be classified according to the magnitude of the SMNS 
spin-down time tgd (expected to vary in the range of weeks to years). For ex- 
ample. X-ray features associated with "light" elements might arise when tgd ~ 
weeks, whereas Fe features could be found when t^d ~ months. 

Simple preliminary estimates concerning the physical viability of this pic- 
ture arc encouraging, but more detailed studies of the formation and evolution 
of supramassive neutron stars and of the pulsar-wind bubbles that they drive 
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are needed. Several unique observational predictions could serve to discriminate 
between this model and alternative scenarios. 

References 

Ballantyne, D.R., et al.: 2002, A&A, 389, L74 

Baumgarte, T.W., Shapiro, S.L.,& Shibata, M.: 2000, ApJ, 528, L29 

Bottcher, M., & Fryer, C.L.: 2001, ApJ, 547, 338 

Chevalier, R.A., & Luo, D.: 1994, ApJ, 421, 225 

Cook, G., Shapiro, S.L., & Teukolsky, S. A.: 1994, ApJ, 424, 823 

Dermer, CD.: 2002, preprint (astro-ph/0211300) 

Esin, A.A., k Blandford, R.: 2000, ApJ, 534, L151 

Fryer, C.L., k Heger, A.: 2000, ApJ, 541, 1033 

Fryer, C.L., k Kalogera, V.: 2001, ApJ, 554, 548 

Gardiner, T.A., k Frank, A.: 2001, ApJ, 557, 250 

Granot, J., & Guetta, D.: 2002, preprint (astro-ph/0211136) 

Guetta, D., k Granot, J.: 2003, MNRAS, in press (GG03) 

Heger, A., Langer, N., k Woosley, S.E.: 2000, 528, 368 

Heger, A., k Woosley, S.E.: 2002, preprint (astro-ph/0206005) 

Inoue, S., Guetta, D., k Pacini, F. 2003, ApJ, 583, 379 

Jun, B.-L: 1998, ApJ, 499, 282 

Kallman, T.R., Meszaros, R, k Rees, M.J.: 2001, preprint (astro-ph/0110654) 
Konigl, A., k Granot, J.: 2002, ApJ, 574, 134 (KG02) 
Livio, M., k Pringle, J.E.: 1998, ApJ, 505, 339 

McLaughlin, G.C., Wijers, R.A.M., k Brown, G.E.: 2002, ApJ, 567, 454 
Meszaros, P., k Rees, M.J.: 1997, ApJ, 482, L29 
Panaitescu, A., k Kumar, P.: 2002, ApJ, 571, 779 
Piro, L., et al.: 2000, Science, 290, 955 

Popham, R., Woosley, S.E., k Fryer, C: 1999, ApJ, 518, 356 

Reeves, J.N., et al.: 2002, Nature, 416, 512 

Ruffert, M., k Janka, H.-Th.: 1999, A&A, 344, 573 

Salgado, M., et al.: 1994, A&A, 291, 155 

Shapiro, S.L.: 2000, ApJ, 544, 397 

Shapiro, S.L., k Shibata, M.: 2002, ApJ, 577, 904 

Shibata, M., Baumgarte, T.W., k Shapiro, S.L.: 2000, PhRvD, 61, 044012 
Spruit, H.C.: 2002, AkA, 381, 923 

van Putten, M.H.P.M., k Ostriker, E.G.: 2001, ApJ, 552, L21 
Vietri, M., k Stella, L.: 1998, ApJ, 507, L45 
Waxman, E., k Draine, B.: 2000, ApJ, 537, 796 

Woosley, S.E., Zhang, W., k Heger, A.: 2002, preprint (astro-ph/0211063) 



